NASA/TM— 1999-000000 



A Mixed Multi-Field Finite Element 
Formulation for Thermopiezoelectric 
Composite Shells 


Ho-Jun Lee 

Glenn Research Center, Cleveland, Ohio 

Dimitris A. Saravanos 

Ohio Aerospace Institute, Cleveland, Ohio 


Prepared for the 

Sixth Annual International Symposium on Smart Structures and Materials 
sponsored by the International Society for Optical Engineering 
Newport Beach, California, March 1-5, 1999 


National Aeronautics and 
Space Administration 


Glenn Research Center 


May 1999 



Available from 


NASA Center for Aerospace Information 
7121 Standard Drive 
Hanover, MD 21076 
Price Code: A03 


National Technical Information Service 
5285 Port Royal Road 
Springfield, VA 22100 
Price Code: A03 



A MIXED MULTI-FIELD FINITE ELEMENT FORMULATION 
FOR THERMOPIEZOELECTRIC COMPOSITE SHELLS 


Ho-Jun Lee 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 

and 

Dimitris A. Saravanos* 

Ohio Aerospace Institute 
Brook Park, Ohio 44142 


ABSTRACT 


Analytical formulations are presented which account for the coupled mechanical, electrical, and 
thermal response of piezoelectric composite shell structures. A new mixed multi-field laminate 
theory is developed which combines “single layer” assumptions for the displacements along 
with layerwise fields for the electric potential and temperature. This laminate theory is 
formulated using curvilinear coordinates and is based on the principles of linear 
thermopiezoelectricity. The mechanics have the inherent capability to explicitly model both the 
active and sensory responses of piezoelectric composite shells in thermal environments. Finite 
element equations are derived and implemented for an eight-noded shell element. Numerical 
studies are conducted to investigate both the sensory and active responses of piezoelectric 
composite shell structures subjected to thermal loads. Results for a cantilevered plate with an 
attached piezoelectric layer are compared with corresponding results from a commercial finite 
element code and a previously developed program. Additional studies are conducted on a 
cylindrical shell with an attached piezoelectric layer to demonstrate capabilities to achieve 
thermal shape control on curved piezoelectric structures. 

INTRODUCTION 


The potential performance advantages from utilizing smart structural components in advanced 
aerospace applications has lead to the extensive recent development of both analytical and 
experimental methods to characterize the behavior of these materials. Crawley (1994), Rao and 
Sunar (1994), and Rogers (1993) provide detailed overviews of smart structure technologies. 
Piezoelectric materials represent one of the more common materials currently being investigated 
for smart structures applications due to their inherent capability to be utilized as both actuators 
and sensors. Of practical interest are curved piezoelectric composite shells, since these 
configurations are most commonly used in various aerospace structural applications. In 
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addition, most aerospace applications involve operations in changing thermal environments. 
Thus, in order to effectively and accurately characterize the response of smart structures with 
piezoelectric materials, the coupling effects from both the curvilinear geometry and 
thermopiezoelectricity must be accounted for in the analysis. 

Numerous theories and models have been proposed to analyze piezoelectric composite beams, 
plates, shells, and solids. Saravanos and Heyliger (1998) provide a comprehensive review and 
classification of these analytical models, consequently only a brief review concerning 
piezoelectric shells is provided herein. Docmeci (1990) investigated the vibration of single 
layered piezoelectric shells. Lammering (1991) developed a Reissner-Mindlin shear deformable 
shell finite element, while Tzou and Gadre (1989) derived equations for thin piezoelectric shells 
based on the Kirchoff-Love hypothesis. Shell formulations and finite elements which explicitly 
account for the piezoelectric coupling have also been reported. Tzou and Ye (1996) modeled 
piezoelectric shells as a layerwise assembly of triangular elements. Heyliger et al. (1 996) 
developed a discrete-layer laminated piezoelectric shell element. Saravanos ( 1 997) implemented 
a mixed shear-layerwise theory for piezoelectric shells. 

Limited research has been performed to assess the implications of thermal effects on 
piezoelectric composite shells. Tauchert (1992) constructed a piezoelectric plate theory for 
thermal behavior based on Kirchoff s assumptions in which the charge equation was not 
considered. Rao and Sunar (1993) developed a general finite element formulation for 
thermopiezoelectric plates. Tzou and Howard (1994) and Tzou and Bao (1995) postulated 
consistent, yet uncoupled, formulations for thin, shallow thermopiezoelectric shells based on 
classical laminate theory assumptions. Chandrashekhara and Kolli (1995) presented an 
uncoupled first-order shear finite element formulation for shallow piezoelectric composite shells 
under thermal loads. Tzou and Ye (1994) developed a three-dimensional finite element for thin 
thermopiezoelectric solids. Tzou and Bao (1997) included large-rotation geometric nonlinearity 
into their previously developed uncoupled models for thermal piezoelectric shells. Lee and 
Saravanos ( 1 996, 1 997) developed electromechanically coupled mechanics and finite elements 
for thermal laminated composite piezoelectric beams and plates, respectively, which use 
layerwise approximations for the displacements, electric potentials and temperature fields. Exact 
solutions for specific configurations of thermopiezoelectric plates and cylindrical shells have 
been reported by Dube et al. (1996a, 1996b), Xu and Noor (1996), Xu et al. (1997) and Choi et 
al. (1997). 

Generally, most of the works described use approximate uncoupled models to account for the 
behavior of thermopiezoelectric laminated shells. However, neglecting the piezoelectric 
coupling that arises from the pyroelectric and thermal expansion coefficients can lead to 
substantial errors in predicting the structural response. Thus, this paper presents a coupled 
piezoelectric shell theory which incorporates thermal effects and has the capability to accurately 
predict the mechanical, electrical, and thermal response of thin and intermediately thick 
piezoelectric composite shells. The coupled behavior is captured at the material level through 
the thermopiezoelectric constitutive equations and the governing equations are formulated and 
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solved in curvilinear coordinates. The mechanics incorporate different types of through-the- 
thickness kinematic approximations for all field variables. First-order shear theory assumptions 
are used for the displacements, while discrete-layer theory approximations are implemented for 
the electric potential and temperature. The combination of mixed through-the-thickness 
approximations for the displacement, electric potential and temperature is a unique feature of the 
"mixed piezothermoelastic shell theory" (MPST) which enables the accurate analysis of thin and 
moderately thick piezoelectric shells. 

THERMOPIEZOELECTRIC LAMINATED SHELLS 

This section describes the analytical formulation for curvilinear thermopiezoelectric laminates 
with embedded sensory and active piezoelectric layers. The curvilinear laminate configuration is 
shown schematically in Figure 1 . Each ply of the laminate remains parallel to a reference 
curvilinear surface A 0 . An orthogonal curvilinear coordinate system O^ is defined, such that the 
axes t, and r| lie on the curvilinear reference surface A 0 , while the axis t, remains straight and 
perpendicular to the layers of the laminate. A global Cartesian coordinate system O xyz is used to 
define A 0 , hence, a point r=(x,y,z) on the curvilinear laminate is, 

*5,tl,0=r o (5,ilK£ (!) 

where, r 0 =(x 0 ,y 0 ,z 0 ) are the Cartesian coordinates of the reference surface A 0 , and C indicates the 
unit vector perpendicular to the reference surface. 



Figure 1 : Curvilinear piezoelectric laminate and coordinate systems. 
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Governing Material Equations 

Each ply is assumed to consist of a linear piezoelectric material with properties defined in the 
orthogonal curvilinear system O^, with constitutive equations of the following form. 


o, = C/S, - e lt E k - X,e 
D i = * E ?A + pfi 


where: i, j = 1,...,6 and k, 1= 1,...,3; o, and Sj are the mechanical stresses and engineering strains 
in vectorial notation; E k is the electric field vector; D, is the electric displacement vector; C,j is 
the stiffness tensor; e /; is the piezoelectric tensor; e tt is the electric permittivity tensor of the 
material; X, is the thermal expansion vector; 0 = AT = T - T c is the temperature difference from 
the thermally stress free reference temperature T 0 ; and p k is the pyroelectric vector. Superscripts 
E and S indicate constant electric field and strain conditions, respectively. The axes 1, 2, and 3 of 
the material are parallel to the curvilinear axes £, r\, and <^, respectively. The materials are 
assumed to be monoclinic class 2 crystals with a diad axis parallel to the z axis (Nye, 1964). The 
assumed material class is general enough, such that Eqs. (2) may encompass the behavior of off- 
axis homogenized piezoelectric plies, as well as the passive composite plies. The tensorial strain 
Sy and electric field components in a curvilinear coordinate system are related to the 
displacements and electric potential, respectively (Saravanos, 1997). 

Mixed Multi-Field Laminate Theory 

The proposed thermopiezoelectric shell theory combines linear displacement fields through the 
thickness of the laminate for the displacements u and v (along the £ and r| axes respectively) with 
layerwise electric potential and temperature fields through the laminate, consisting of N discrete 
continuous segments (see Figure 2). Previous works by Lee and Saravanos (1996, 1997) on 
layerwise theories for thermal piezoelectric beams and plates have demonstrated the advantages 
and necessity of layerwise approaches for capturing the complex electric fields and interactions 
which are present in piezoelectric actuators and sensors. Consequently, the present mixed 
piezothermoelastic shell theory (MPST) will entail the capability to (1) accurately and efficiently 
model thin and moderately thick laminated piezoelectric shells with arbitrary laminations and 
electric configurations, (2) capture the through-the-thickness electric heterogeneity induced by 
the embedded piezoelectric sensors and actuators, and (3) represent the changing thermal 
gradients through the laminate plies. 

The discrete-layer kinematic assumption effectively subdivides the laminate into N-l 
sublaminates (or discrete-layers). A continuous linear electric potential and temperature variation 
is assumed in each sublaminate, such that a C„ continuity results through-the-thickness of the 
laminate (see Figure 2). In this formulation, the subdivision can be arbitrarily altered according 
to the configuration of the piezoelectric layers or the desired level of approximation. Thus, the 
resulting “zig-zag” shape of the approximation can be a-priori controlled to vary the detail of 
approximation from a single-layer to a finely graded field. The various field variables are 
approximated by the following form, 
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«&n.C.O= w°(£,i 1,0 + CP^tlO 

v(£,r|,C,0= v°(^,T|,0 + CP^ry) 

w(&,r\,0)= w o (§,ii,0 

<N 5 .ti,C, 0 = E 4 /( 5 ,vWO 

J l 

0 (§,n.C,o= E e'&rymo 

y=l 


where u°, v°, w° are displacements along the p and C, axes, respectively, on the reference 
surface A 0 ; superscript] indicates the points Q at the beginning and end of each discrete layer; (J) 1 
and 9' are the electric potential and temperature at each point Q, respectively (see Figure 2); 'P(z) 
are interpolation functions; and (3, and P n are the rotation angles defined as. 




O 

Si i 





O 

S 22 



( 4 ) 


where R; are the local radii of curvature (Figure 1). Linear interpolation functions 'F(z) are 
considered in this paper. 



(a) Smart laminate 


(b) Piezoelectric shell laminate theory assumptions 


Figure 2: Typical piezoelectric laminate configuration, 
(a) Concept; (b) Assumed through-the-thickness fields. 
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The Love's assumption for shallow shells is further implemented, i.e. the local radii of the shell 
are substantially higher than the thickness (h/R, « 1 ), yielding ( 1 + C/Rj = 1 ). In the context of 
Eqs. ( 3 ), the engineering strains become 

S'tt, J\,U) = S°a,T),t) + Ckfc,r\ ,t) /= 1 , 2,6 

S 2 (£,,r\,(,t)= 0 (5) 

Sfc, ri,C,0= S°& ry) i=4,5 


where S° and k are the strain and curvatures at the reference surface, defined as follows: 


iO 1, o Si |,ij 0 - W 0 

— 1 1 1 -±- 'll 


-« + — i v °) 

O * () /? 

£| 1 &22 K 


O o I/O Si 2.£ W 0 

= — (v n + — " > + y 

«22 i’ll ! 


O O I/O O l 1 ,T) (ji i/o ^i.2,,1. /ii 

5 6 = — ( v £- ) + — (a,-— V°) 
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( 6 ) 
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( 7 ) 


where R, are the local radii of curvature, and g°, ,, g ° 22 are the components of the metric tensor on 

the surface A 0 (£= 0 ) defined as g° = /x„V>i +z i W £22 • 

The electric field vector also becomes, 

£,.(£, T1,C,0= E £,U,Tl,0^(O i-1,2 


7 = 1 
JV 


£ 3 (S,r|,C, 0 = E 

7 = 1 


( 8 ) 


where {E J } is the generalized electric field vector defined as: 
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( 9 ) 


£/=- — i Ei = -^- E{=-& 

o o 

£ll S 22 


Equations of Motion 

The variational form of the equations of motion in the orthogonal curvilinear system is, 

| { -bH(S,E)-bT }\J\(%dr\dZ + J bu i T x i + JbtyDdS = 0 
' r T r „ 


( 10 ) 


where 8H is the variation of the electric enthalphy and 8T the variation of the kinetic energy, 
defined as. 


8 H = 5S.o-6E.Di 

i i k k 

8 T - -5up j u j 


5SC S -5Se ErbSXQ-bE.e.S -bE.z.E.-bE'pfi 

i ij j 1 1/ l t 1 k kj j k kl l k 


(ii) 


where xand D are the surface tractions and electric displacement acting on surfaces T s and r D , 
respectively 

V<YV ® =D j”j .3 (12) 

[J] is the Jacobian matrix of the transformation between the global cartesian and curvilinear 
system. For the curvilinear system the determinant of the Jacobian takes the form: 

MH ^nO l = = (13) 


It is now possible to rearrange the variational statement and separate the through-the-thickness 
integration, as follows: 


~(bH L (S,E) - 8 T l ) d\dx\ + £ (bu t T x ^D)dT=0 


(14) 


where A c is the curvilinear reference surface, x and D are the surface tractions and electric 
displacement on the boundary surface T; 8H L and 5T L are the variations of the electric enthalpy 
and kinetic energy of the laminate, defined as 

<bH L ,bT L > = f h <bH,bT>g x \g° 2 d( (15) 

0 
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and h is the laminate thickness. By combining Eqs. (11) and (15) and integrating through-the- 
thickness, the variation of the electric enthalpy of the piezoelectric laminate is obtained as a 
quadratic expression of the generalized strain, electric field, and the generalized laminate 
matrices. 


8 H l - bS°A IJ S°>?,S°B, J k 1 ^k ] B.S°+bkp iJ k J - 

~t (,bS;E2E^&E^E2spUE~E, m *6E’E’k) 'EE bE’G^E," 

1 m =l n = 1 I 16 ) 

( 6S°A "6" + 60 m A '"S,° + 6^A“0 m +60 m A / m A.) E ^>E k m P" n d m 


N N 


m - 1 


m = 1 * = I 


In the above equation, [A], [B], and [D] are the stiffness matrices; [E m ] overbar and overhat are 
the piezoelectric matrices; and [G mn ] are the matrices of electric permittivity of the curvilinear 
laminate described by Saravanos (1997). 


The generalized thermal matrices are new, and include the thermal expansion matrices 
[A m ] overbar and overhat. 


< A"', A'” > 


suft°2 E f c ‘" vr"(0<i,wc, 

m h, 


M, 2 


(17) 


and the pyroelectric laminate matrices [P mn ], 

C - E f ( " [ />, 'r"(o , F”(f)rfc , k-- 1,2 

h\ J C , 

- *il 8 a E f ^^"(CmOdC 

/=! J C, 


(18) 


where L is the number of plies in the laminate. 


The kinetic energy of the laminate takes the form, 


/'=!,. ..,3, y=l,2 


(19) 


where u°— {u°,v°,w 0 } and P-{P^, P n }; p A , p B , p D are the generalized densities, expressing the 
mass, mass coupling and rotational inertia per unit area of the laminate (Saravanos, 1997). 

Finite Element Formulation 

The previous formulation of the governing equations in the orthogonal curvilinear system and the 
generalized variational statement in Eq. (11) enables the development of structural solutions by 
using approximations for the generalized multi-field state variables (displacements, rotation 
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angles, electric potential and temperature) of the following type, 

M 

E «>"(') tf'&T]), 7 = 1,..., 3 

M 

M 

p,&ty)= EPytO^^Ti), 7=1,2 

* ( 2 °) 

<|>"(5,T17)= Er'W^'^Tl), m = l,...,^v 

M 

M 

e m (£,ri,o= £ e m '(/)^^,Ti), «=i,...jv 

/= i 

where superscript i indicates the reference surface displacement, rotation angle and generalized 
electric potential components corresponding to the i-th in-plane interpolation function N'(£,,r|). 
Implementation of local interpolation functions in Eqs. (19) results in the finite element based 
solutions. Combining Eqs. (6), (7), (9), and (20) results in the strain, electric field and 
temperature interpolation matrices. Substituting into the variational equations of motion Eqs. 
(14), and collecting the coefficients as mandated by Eqs. (16) and (19), .the governing dynamic 
equations of the structure are expressed in the discrete matrix form, 

\mj oil <«)1 J[*j [Cl]| Ml j MOHOVHK.eim 1 

Submatrices K uu , K u<1 and K w indicate the elastic, piezoelectric and permittivity matrices; K u0 and 
are the thermal expansion and pyroelectric matrices of the structure; and M uu is the mass 
matrix. Superscripts F and A indicate the partitioned submatrices in accordance with the sensory 
(free) and active (applied) electric potential components, respectively, as explained previously 
(Lee and Saravanos, 1997). Thus, the left-hand side includes the unknown electromechanical 
response of the structure {u, <j> F }, i.e. the resultant displacements and voltage at the sensors. The 
right-hand includes the excitation of the structure in terms of mechanical loads, applied voltages 
(J> A on the actuators, applied temperature loads and electric charge at the sensors Q F (t). 

Based on the above formulation, an eight-noded (M=8) finite element was developed with bi- 
quadratic shape functions of the serendipity family. Reduced integration was applied selectively 
on the calculation of shear elastic and piezoelectric stiffness terms, as this has been found to 
improve the accuracy of the numerical predictions. These discrete equations can be easily 
condensated into the following independent equations for the sensory electric potentials, 

{<t> f f = [<£] 1 ([*«] d)"! * I**,] {0} -i*Ti m - ie F (')i ) (22) 

and the structural displacements. 
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( 23 ) 


mm - < kJi k£i 1 [<J] * iO ) (V> + 
tOO' 1 iB F m -aO[<*]‘'[^e] *[^ e ])(0} 

NUMERICAL RESULTS AND DISCUSSION 


Evaluations of the developed mechanics and applications to various composite structures with 
piezoelectric actuators and sensors are presented in this section. In order to represent the location 
of the piezoelectric layers, the standard laminate notation is expanded, such that piezoelectric 
layers are indicated by the letter p. The properties of both the piezoelectric and composite 
material are listed in Table 1. The structural configurations which are studied include a clamped 
plate and a circular cylindrical shell. 

Clamped Plate 

The first case examined consists of a 50.8 mm x 25.4 mm x 5.08mm [0 8 /p] plate composed of 
carbon/epoxy with an attached piezoceramic layer as shown in Figure 3. The plate is clamped on 
one side and is subjected to a uniform thermal load of 50°C. The objective of this part of the 
numerical study is to verify the accuracy of the current formulation for both the sensory and 
active response with corresponding results from a 20-node quadratic piezoelectric continuum 
element (MSC/ABAQUS, 1996), as well as a previously developed 4-node plate element (Lee 
and Saravanos, 1996) that implements a full layerwise laminate theory. To the author’s best 
knowledge, the continuum element neglects pyroelectric effects in the piezoelectric formulation, 
so the pyroelectric constants were intentionally neglected in both the shell and plate program 
analysis during the comparisons with the continuum element. Subsequent studies are performed 
using only the plate and shell elements to investigate the influence of incorporating the 
pyroelectric effect. Since the three programs utilize different formulations and elements, each 
analysis required different finite element meshes. The shell program utilized a 1 0x5 mesh with 
two discrete layers (one for the carbon/epoxy and one for the piezoelectric). A finer 20x10 mesh 
was used for the 4-noded plate element with two discrete layers. The continuum element 
analysis incorporated a 10x5x2 mesh, with the two elements through the thickness of the plate 
corresponding to the two discrete layers used in the shell and plate programs. All displacements 
(w) and electric potentials (4>) presented in this section are along the centerline of the plate 
(y/b=0.5) and are nondimensionalized using the following relationships: 

w ' = lOOw / h ; ({)' = 4> d 3| 10 5 / h (24) 

where h represents the plate thickness and d 3 , represents the piezoelectric charge constant. 
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TABLE 1: MATERIAL PROPERTIES 



Piezoceramic 

Carbon/Epoxy 

Elastic Moduli (10 9 Pa): 

E u 

69.0 

142.0 

E 2 2 

69.0 

10.3 

E33 

69.0 

10.3 

^12 

26.5 

7.20 

G b 

26.5 

7.20 

®23 

26.5 

4.29 

Poisson’s Ratio: 

l) ]2 

0.30 

0.27 

U 13 

0.30 

0.27 

U 23 

0.30 

0.20 

Thermal Expansion Coefficient ( 10' 6 / °C): 

“n 

1.20 

-0.90 

tt 22 

1.20 

27.0 

Piezoelectric Charge Constant (10 12 m/V): 

^31 

-154.0 

0.0 

Electric Permittivity (10‘ 9 farad/m): 

6 11 6 22 e 33 

15.05 

0.0266 

Pyroelectric Constant (10‘ 3 coulomb/m 2 °C): 

P3 

-0.2 

0.0 

Reference Temperature (°C): 

T 

a 0 

20.0 

20.0 
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Figure 3: Geometry of a clamped plate with an attached piezoelectric layer. 


Sensory Response. The response of the plate with the piezoelectric layer operating as a sensor 
is presented in this section. In the sensory mode, the piezoelectric layer is free to develop a 
corresponding electric potential in response to the thermally induced deformation of the plate 
through the direct piezoelectric effect. The resulting displacements predicted by the three 
different codes show almost exact agreement as depicted in Figure 4. The corresponding sensory 
electric potential which develops in the piezoelectric layer is shown in Figure 5. Although the 
continuum element electric potential results tend to be slightly higher due to the consideration of 
the through-the-thickness Poisson’s effect, there is good overall agreement between the three 
elements. 

Active Response. The piezoelectric layer can also be used as an actuator by applying a voltage 
differential to the piezoelectric layer and utilizing the converse piezoelectric effect to alter the 
thermal distortion of the plate. Figure 6 shows the centerline deflection of the plate for three 
different applied electric potentials (4> A *). The three codes show almost exact agreement and 
demonstrate the capability to achieve thermal shape control by applying electric voltages. 
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x/a 


Figure 4: Thermally induced sensory displacements of a [0 g /p] clamped plate. 



Figure 5: Thermally induced sensory electric potentials of a [0 8 /p] clamped plate. 
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Figure 6: Displacements under active electric potentials of a [0 8 /p] clamped plate. 


Pyroelectric Effects. The effect of incorporating the pyroelectric constant (p 3 ) into the analysis 
is presented in this section. The previous two cases have neglected the pyroelectric effect in 
order to compare with the continuum element. Results for the shell and plate elements are 
presented with the piezoelectric layer in a sensory mode. Figures 7 and 8 show the centerline 
displacements and the centerline electric potential, respectively. The case with p 3 =0 represents 
the case without considering pyroelectric effects shown in Figures 4 and 5. Both figures show 
the almost exact agreement between the plate and shell elements, as well as the significant 
impact that pyroelectric effects can have on the displacement and electric potential of a smart 
structure operating in changing thermal environments. This significant influence on the 
structural response is caused by the coupling introduced by the piezoelectric layer and can be 
modeled only through coupled formulations, such as the present one. Uncoupled formulations 
are unable to capture this effect and can result in significant errors when modeling thermal 
piezoelectric structures. 
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Figure 7: Pyroelectric effects on sensory displacements of a [0 g /p] clamped plate. 
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Figure 8: Pyroelectric effects on sensory electric potentials a [0g/p] clamped plate. 
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Cylindrical Shell 


The second problem studied consists of a [0 9 /p] (where 0° is along the £ axis) circular cylinder 
with a radius R=0.76 m, a L=1 .524 m, and a thickness h=7.62 mm as depicted in Figure 9. Due 
to the symmetry of the problem, only Ve of the cylinder is modeled using a 10x10 mesh with two 
discrete layers (one for the carbon/epoxy and one for the piezoceramic). A clamped-clamped 
configuration in which both ends of the cylinder (z/L = 0 and z/L = 1) are fixed is examined. The 
cylinder is subjected to three different types of thermal loads (a uniform thermal load and two 
types of sinusoidally varying temperatures along the hoop direction) as shown in Figure 10. 
Results from the current shell element are presented to investigate the influence of the curved 
structure on achieving thermal shape control under the different boundary conditions and applied 
thermal loads. All displacements (w) and electric potentials ((j>) presented in this section are 
along the hoop direction at z/L = 0.5 and r/R =1 .0 of the cylinder and are nondimensionalized 
using the following relationships: 

w ' = 1 00 w / h ; (j) ’ = (J) d 3 , 1 0 4 / h (25) 

where h represents the cylinder thickness and d 3l represents the piezoelectric charge constant. 


x 



Figure 9: Geometry of a [0 g /p] circular cylinder with an attached piezoelectric layer. 
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T - 50 °C 


T = 50 °C 


T = 25 °C 


(a) Uniform T (b) T cos y ( c ) T cos 2y 


Figure 10: Different types of applied thermal loads. 

Sensory Response. The sensory response of the cylinder for the three different applied thermal 
loads (Figure 10) are presented in Figures 11-13. Each figure contains results for the sensory 
displacement (w*) and electric potential (4)*) for two cases to show the effects of incorporating 
the pyroelectric constant in the analysis (the p 3 =0 case neglects pyroelectric effects). Figure 1 1 
shows the displacement and sensory electric potential induced from applying a uniform thermal 
load of 50°C. A uniform deflection of the cylinder is achieved, which maintains the original 
circular shape. The sensory electric potential also shows a uniform response that corresponds to 
the displacement. The incorporation of pyroelectric effects produces larger displacements and 
electric potentials. The results from applying a cosine varying temperature are shown in Figure 
12. The deflection of the cylinder now shows a sinusoidal pattern that translates into either an 
oval-shaped or figure eight shaped circumferential displacement, depending on whether 
pyroelectric effects are modeled. The electric potential only displays a sinusoidal pattern when 
pyroelectric effects are included and becomes almost zero when the pyroelectric effect is 
neglected. Figure 13 shows the results of applying a double cosine varying thermal load. A 
sinusoidally varying displacement and electric potential pattern is produced which corresponds to 
an oval-shaped deflection of the cylinder. Once again, incorporating the pyroelectric effect leads 
to increased displacements and electric potentials. 
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Figure 11: Sensory response of a [0,/p] circular cylinder under a uniform thermal load. 
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Figure 13: Sensory response of a [(yp] circular cylinder under a double cosine thermal load. 


Active Response. The active response of the cylinder under the three applied thermal loads is 
examined in Figures 14-16. Each figure shows the displacement of the cylinder for three cases: 
the sensory thermal deflection resulting when the electric potential of the piezoelectric layer 
remains free, a grounded thermal configuration (zero electric potential is applied on both 
piezoelectric terminals), and an active thermal deflection when a non-zero electric potential is 
applied. The magnitude of the applied electric potential remains the same for all three cases, 
although the form of the electric potential varies to correspond to the type of thermal load. 

Figure 14 shows the deflection of the cylinder under a uniform thermal load (50 °C). Application 
of increasing uniform applied electric potentials produces a noticeable decrease in the sensory 
thermal deflection. Figure 15 depicts the displacement of the cylinder under a cosine varying 
thermal load and shows that although some changes in the overall shape of the cylinder can be 
obtained by applying active electric potentials, no significant overall reduction of the thermal 
deflection is achieved. Figure 1 6 illustrates the deflections produced under a double cosine 
variation of the thermal load. For this case, the double cosine applied electric potential 
effectively minimizes the thermal deflection. 
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Figures 14-16 also indicate that a significant reduction in thermal deflection occurs between the 
sensory and grounded ((j) A = 0) configurations. This reduction is comparable to the thermal 
deflection achieved by applying the non-zero electric potentials and is attributed to the influence 
of pyroelectric effects on the response of the cylinder. The modeling of this effect is a direct 
consequence of the coupled formulation. Under thermal loads, the pyroelectric effect has a 
significant impact on the sensory electric potential, cJ> F , in Eq. (22). Incorporating this increased 
sensory electric potential into the structural displacements, Eq. (23), produces a corresponding 
effect on the deflections. However, when the piezoelectric material is in an active configuration, 
the sensory electric potentials are eliminated, significantly reducing the pyroelectric effect on the 
displacements as shown in Figures 14-16. While these results help illustrate the advantages of 
the current coupled formulation, it also shows a potential drawback of using piezoceramics as 
sensors in high temperature applications. Although typical piezoceramics have higher use 
temperatures than piezopolymers, the pyroelectric effect can lead to significant changes in the 
mechanical state of the structure and appropriate design considerations will be required for 
effective operations in thermal environments. 
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Figure 14: Displacements of a [CVp] circular cylinder under a uniform thermal load. 
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Figure 15: Displacements of a [Og/p] circular cylinder under a cosine thermal load. 
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Figure 16: Displacements of a [0 9 /p] circular cylinder under a double cosine thermal load. 
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SUMMARY 


A multi-field laminate theory and mechanics to model the coupled mechanical, electrical, and 
thermal behavior of piezoelectric composite shells were described. The corresponding equations 
of motion were developed in curvilinear coordinates and implemented into a quadratic shell 
element. Analytical studies were conducted to compare the accuracy of the current formulation 
with a commercial finite element code, as well as with a previously developed formulation. 
Additional numerical studies have shown the significant impact of piezoelectric coupling and 
pyroelectric effects on the sensory response of thermal piezoelectric structures. The response of 
cylindrical shells under combinations of applied thermal and electric loadings was quantified and 
numerical examples illustrate the feasibility of active thermal shape control on curved structures. 
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